Disorder-Induced Shift of Condensation Temperature for Dilute Trapped Bose Gases 
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We determine the leading shift of the Bose-Einstein condensation temperature for an ultracold 
dilute atomic gas in a harmonic trap due to weak disorder by treating both a Gaussian and a 
Lorentzian spatial correlation for the quenched disorder potential. Increasing the correlation length 
from values much smaller than the geometric mean of the trap scale and the mean particle distance 
to much larger values leads hrst to an increase of the positive shift to a maximum at this critical 
length scale and then to a decrease. 
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The experimental realization of Bose-Einstein conden- 
sates (BECs) in magnetically trapped dilute atomic va- 
pors [1, 2] has renewed the theoretical interest in the 
phase structure of ultracold Bose gases. In particular, 
many research groups have investigated how the transi- 
tion temperature T c changes due to a weak local repul- 
sive two-particle interaction V^ mt ^(x) = Aira/M <5(x) in 
different trap configurations. Here a denotes the s-wave 
scattering length, M the atomic mass, and we employ 
units with fee = h = 1. Theoretically most demanding is 
the homogeneous case [3] where the Bose gas of particle 
density n resides in an infinitely wide trap with a flat 
potential. There the leading shift AT C /T C (0) = dan 1 / 3 
with respect to the interaction-free critical temperature 
T c (0) = 27r[n/C(3/2)] 2/3 /M is completely due to long- 
wavelength, critical fluctuations that have to be de- 
scribed nonperturbatively. It can be shown that this 
shift follows from the classical limit of many-body the- 
ory, where only the zero Matsubara modes of the Bose 
fields are included [4], A self-consistent variational re- 
summation [5] of the resulting infrared divergent series 
for c\ yields the value 1.27 ± 0.11 [6] which agrees with 
the Monte Carlo calculations 1.32 ± 0.02 [7, 8]. In con- 
trast to that, harmonic trapping potentials remove the 
critical long- wavelength fluctuations and reduce the frac- 
tion of atoms taking part in nonperturbative physics 
at the transition point. As a result, the leading shift 
AT C /Tc°^ = dia/X^ in the critical temperature with 
respect to its noninteracting value — [7V(jj g /£(3)] 1//3 
is negative and can be calculated by simple perturba- 
tive methods, for instance a Hartree-Fock mean-field ap- 
proximation [9]. Here a£ 0) = \J 2n / [t!; 0) M] denotes the 
thermal wavelength and N the boson number, whereas 
uj g = (uJxUJyLdz) 1 / 3 is the geometric mean of the respec- 
tive trap frequencies. The analytically known constant d\ 
has the numerical value d\ = —3.426 which could be con- 
firmed in a recent experiment within the error bars [10]. 
Note that a theoretical investigation of the crossover be- 
tween both physically distinct cases of the homogeneous 
and the trapped Bose gas has been elaborated in Ref. [11] 
with the help of general power-law potentials. Further- 



more, corrections to the leading interaction-induced shift 
in the critical temperature have been investigated in the 
Rcfs. [12-15]. 

Apart from the contact interaction, another interesting 
physical parameter for ultracold atomic vapors is frozen 
disorder. It can be created artificially by laser speck- 
les [16-19], incommensurate lattices [20], or different lo- 
calized atomic species [21, 22]. However, random po- 
tentials can also arise naturally via the spatial fluctu- 
ations of the electric currents generating the magnetic 
wire traps [23, 24]. Thus, one can raise the fundamental 
question how the Bose-Einstein condensation tempera- 
ture is affected by weak disorder. To this end we consider 
the quite commonly employed model of a Gaussian dis- 
tributed disorder potential V(x) with the ensemble aver- 
ages V(x) = and V(x)V(x.') = i?(x— x'). In particular, 
we are interested in correlation functions i?(x— x') which 
decay with a characteristic correlation length £. A typical 
example for such a correlation function is a Lorentzian 
distribution, which occurs in atom chips [25], more com- 
plicated correlation functions can be designed by using 
laser speckles [26] . Note that the latter experiments allow 
to tune the corresponding correlation length £ in a wide 
parameter range by changing, for instance, the distance 
between the scattering surface and the magneto-optical 
trap. So far, the disorder-induced shift of the critical 
temperature has only been analyzed for a homogeneous 
Bose gas in the limit £ — * of a short-range correlation 
.R(x) = R5(x). A perturbative calculation of Lopatin 
and Vinokur [27], motivated by studying the superfluid 
motion of helium in porous media such as Vycor [28-30] , 
yields 
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The negative sign for the homogeneous Bose gas is com- 
patible with the fact that disorder leads to an effective 
attractive interaction when the replica method is applied 
[31]. 

The purpose of this note is to determine the leading 
disorder-induced shift in the critical temperature for a 
weakly interacting Bose gas in an external harmonic 
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trap U(x) with an additional weak disorder potential 
which is characterized by an arbitrary correlation func- 
tion i?(x — x'). According to our introduction we expect 
that disorder leads to an opposite sign of the shift of the 
critical temperature than repulsive interaction. Since dis- 
order appears within the replica method as an effective 
attractive interaction [31], the particle density in the cen- 
ter of a harmonic trap is increased. Thus, there should be 
a positive disorder-induced shift for such a harmonic con- 
finement. Moreover, we should be able to determine it re- 
liably by a perturbative calculation. In first-order pertur- 
bation theory, the effects of both interaction and disorder 
are additive. Therefore, we start with the imaginary-time 
action of a non-interacting Bose gas in both a trap and 
a random potential in standard notation 
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and treat the disorder potential V(x) perturbatively. We 
aim at determining the particle number N — —dil/du 
from the grand-canonical free energy = — lnZ//3, 
where the partition function is determined by the func- 
tional integral 



2?V* f Tty e 



(3) 



For an ideal Bose gas a semiclassical treatment of the 
trap potential leads to 



(4) 



where ( a (z) — J2^=i z n /n a is the polylogarithmic func- 
tion. Correspondingly the first-order correction reads 



= ~h ll dTl ll dT2 J d3xi J Sx2 

Xj R(xi - x 2 )G^(xi, n; x 2 , r 2 )G^(x 2 , r 2 ; xi, n) , (5) 

where the bosonic Green function has the semiclassical 
Fourier-Matsubara decomposition 

G <°>( Xl ,n;x 2 , T2 ) = I t l§f 

m— — oo v ' 
xe -^ m( r 1 -r 2 )+ l k(x 1 -x 2 ) G (0) ^ ro;k; Xl + X2 ^ ( 6 ) 

with G(°)(w ro ,k;x) = [-iuj m + k 2 /2M + U(x)-fi]~\ 
x = (xi+X2)/2, and uj m — 2mrn/ (3. Evaluating in (5) the 
Matsubara series with the Poisson summation formula 
[33] and the Fourier integrals with the formula of Plemelj 



[34], we obtain 
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This disorder correction to the grand-canonical free 
energy is now specialized for a harmonic trap 
t/(x) and for typical correlation functions i?(x) = 
J c? 3 fc/(27r) 3 e lkx i?(k). For instance, in case of a Gaus- 
sian i?c(k) = i?e~ k2 £ 2 / 2 and a Lorentzian RlQs) = 
i?/(l + k 2 ^ 2 ), we find 
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where we have used for the latter the probability integral 
$(z) = 2/0F/ o ^e-* 2 . 

As follows from the theory of critical phenomena [35, 36], 
the phase transition to a Bose-Einstein condensate occurs 
when the spatio-temporal integral over the correlation 
function 

G(xi,ti;x 2 ,t 2 ) 

= <j)Vip* <j)Viptl>(x 1 ,T 1 )^*(x 2 ,T 2 )e- A &*^ (10) 

diverges, i.e. G(w m ,k;x) becomes infinite for oj m = 0, 
k = 0, and some x. This criticality criterion is equivalent 
to the vanishing of its inverse 

G _1 (0,0;x) = G^ _1 (0,0;x) - E(0,0;x) , (11) 

where the first-order perturbative contribution to the 
self-energy is given by 

S( 1 )(x 1 ,T 1 ;x 2 ,r 2 )= J R(x 1 -x 2 )GW(x 1 ,r 1 ;x 2 ,r 2 ) .(12) 

Therefore, we define the critical chemical potential u c 
from the condition 

^ n G- 1 (0,0;x) = 0. (13) 

Up to first order in the disorder we obtain a result which 
does not depend on the trap frequencies: 
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FIG. 1: Disorder-induced shifts (15) and (16) of the critical 
temperature AT C with respect to the interaction-free critical 
temperature T c for the Gaussian (solid) and the Lorentzian 
(dashed) correlation function together with their asymptotics 
(17) and (18) as a function of the dimensionless correlation 
length. 



We conclude from (14) that Bose-Einstein condensation 
occurs with disorder earlier than in the clean case. Note 
that the critical chemical potential for the onset of Bosc- 
Einstein condensation follows from condition (13) also in 
other physical situations, for instance, when anisotropic 
traps are used for rotating Bose gases [37] or when dipo- 
lar gases are considered where the two-particle interac- 
tion contains both a short-range, isotropic contact po- 
tential and a long-range, anisotropic dipole-dipole inter- 
action [38, 39]. 

Now we combine our previous results and determine the 
disorder-induced shift in the critical temperature. To this 
end we evaluate N — N(fi) in the limit /i | by taking 
into account the Robinson formula [40] for the polyloga- 
rithmic function in (4). For the Gaussian correlation we 
find 
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whereas we obtain for the Lorentzian correlation 
AT,.r, R ^ f 1 [C(3)] 1/6 V 



(15) 



x exp 
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Here we have introduced dimensionless units for the cor- 
relation length £ = £,/l os as well as for the disorder 
strength R = M 2 l os R, where l os = l/-y/Mv e denotes 
the oscillator length. Note that both disorder-induced 



FIG. 2: Normalized shift of the critical temperature 
AT Ci g/(Tc°'-R) of a dilute 85 Rb gas in an anisotropic har- 
monic trap with geometric mean frequency u) g = 2n ■ 40 Hz 
and a Gaussian disorder correlation with correlation length 
£ = 10 /im versus particle number N. 



shifts (15) and (16) have similar asymptotics in the limit 
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Figure 1 shows how the shift results (15) and (16) vary 
with the dimensionless correlation length £. Whereas for 
small £ the oscillator length ^ os dominates, the average 
particle distance l p = los/N 1 / 3 becomes important for 
large £. Both length scales meet at yJl os l v = los/N 1 ^ 6 , 
i.e. t^N 1 / 6 = 1, which is roughly where the maximal shift 
occurs. This result will certainly help to design an exper- 
iment which aims at detecting a disorder-induced shift of 
the critical temperature. 

Furthermore, it is clear from (15), (16), and Figure 1 that 
the shifts for Gaussian or Lorentzian spatial correlations 
are qualitatively and even semi-quantitatively the same. 
This finding justifies to apply our results also to other 
correlation functions, which decay with a characteristic 
correlation length £ as, for instance, the laser speckle 
correlations in the Florence experiment [16, 18]. There 
about N = 3 • 10 5 atoms of the 85 Rb isotope reside in 
an anisotropic harmonic trap with the geometric mean 
frequency uj g — 2ir ■ 40 Hz, so that its oscillator length is 
l os = 1.72 /im. The disorder potential is characterized by 
a length scale £ = 10 /mi, i.e. £ = 5.88, and a strength 
R which is measured in units of a frequency u>. As R has 
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the physical units (energy) 2 (length) 3 , the dimensionless 
disorder strength can be estimated by R = (cj/cj g ) 2 ^ 3 . 
A weak disorder potential, which only slightly perturbcs 
the condensate in the Florence experiment, corresponds 
to the frequency u — 2ir ■ 50 Hz, yielding R = 200. In 
Figure 2 the normalized shift of the critical temperature 
AT C /(T C (0) J R) is plotted as a function of the particle num- 
ber N for a Gaussian correlation, so we conclude that 
the disorder-induced shift in the critical temperature is 
in this case of the order of 3 %. With this it is about the 
same size as the interaction-induced shift [9] and should 
thus be measurable [10]. Furthermore, we read off from 
Figure 2 that this shift increases with a decreasing num- 
ber of particles in the Bose gas. 

In this letter we have determined the disorder-induced 
shift of the critical temperature for a Bose gas in a har- 
monic trap. By treating the analytically solvable cases 
of a Gaussian and a Lorentzian correlation function, we 
have estimated that this disorder-induced shift should 
be measurable in laser speckle experiments. Note that 
applying our general perturbative results (5) and (14) 
to a homogeneous Bose gas and to a vanishing corre- 
lation length £ = reproduces the seminal Lopatin- 
Vinokur shift (1). However, in this case the interaction- 
induced shift of the critical temperature suggests that 
additional nonperturbative contributions to the disorder- 
induced shift will certainly occur. First indications for 
this conjecture are provided by a recent momentum-shell 
renormalization group calculation [32] . 
This work was supported by the SFB/TR 12 of the Ger- 
man Research Foundation (DFG). 
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